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1 .Introduction. 



Since 1952 wlien A.C Turing publislied the remarkable paper |T], which a revolutionary 
idea about mechanism of morphogenesis (the development of structures in an organism during 
the life) has been proposed in, nonlinear reaction-diffusion systems of the form 

>^iUt = u^^x + F{u,v), 
\2Vt = + G{u,v), 

have been extensively studied by means of different mathematical methods, including group- 
theoretical methods (see [21 El S] and the papers cited therein). In system ([T]), F and G 
are arbitrary smooth functions, u = u{t,x) and v = v{t,x) are unknown functions of the 
variables t, x, while the subscript t and x denotes differentiation with respect to this variable. 
Notably nonlinear system ([T|) generalizes many well-known nonlinear second-order models used 
to describe various processes in physics [5], biology [6l [7] and ecology [8J. 

In the present paper, we shall consider the diffusive Lotka-Volterra (DLV) system 

XiUt = Uxx + u{ai + biu + civ) , . . 

A2ft = + V{a2 + b2U + C2V), 

which is the most common particular case of reaction-diffusion (RD) system ([T|). System ([T]) 
is the simplest generalization of the classical Lotka-Volterra system that takes into account 
the diffusion process for interacting species (see terms Uxx and Vxx)- Nevertheless the classical 
Lotka-Volterra system was introduced by A.J.Lotka and V.Volterra more than 80 years ago, its 
natural generalization ([2D is still studied because this is one of the most important mathematical 
models. Lie symmetries of have been completely described in ^ (note those can be extracted 
from more general results presented in [21 [3]). 

The problem of construction of Q-conditional symmetries (non-classical symmetries) for ([1]) 
is still not solved even in the case of DLV system ([2|). Moreover, to our best knowledge, there are 
only a few papers devoted to the search of conditional symmetries of systems of PDEs. Notably, 
some general results about Q-conditional symmetries of RD systems with power diffusivities of 
the form 



have been obtained in the recent paper [TU]. However, the results obtained in [TU] cannot be 
adopted for any system of the form ([Tj) because the case / = A; = is a very special and wasn't 
studied therein. 

It should be noted that there are many papers devoted to the construction of such symme- 
tries for the sca/ar non-linear reaction-diffusion (RD) equations of the form [Tn[T2l[T3l[Til[T5l[T6] 



and ([ID with the convective term B{U)Ux (here B{U), D{U) and F{U) are arbitrary smooth 
functions) [HI IH [19] . 



ut = {u^Ux)x + F{u,v) 
vt = {v^Vx)x + G{u,v) 



(3) 



U, = [D{U)U,l + F{U) 



(4) 
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It is well-known that conditional symmetries can be applied for finding exact solutions of 
the relevant equations, which are not obtainable by the classical Lie method. Moreover the 
solutions obtained in such a way may have a physical or biological interpretation (see, e.g., 
examples in fl8\ IT9| 1201 [21] ) what is of fundamental importance. 

The paper is organized as follows. In Section 2, we present two definitions of Q-conditional 
symmetry in the case of RD system ([2]) and show how they are connected with non-classical 
symmetry. In Section 3, we present a complete description of Q- conditional symmetries of the 
DLV system (|2]), i.e. the system of the determining equations for constructing Q-conditional 
symmetries of system ([1]) is derived and analyzed. Here the main theorems presenting these 
symmetries in explicit form are proved. In Section 4, the Q- conditional symmetries obtained 
are applied to reduce the corresponding DLV systems to the systems of ordinary differential 
equations (ODEs) and constructing exact solutions. The properties of an exact solution are 
examined with the aim to provide the relevant interpretation for population dynamics. 

Finally, we present some conclusions. 

2. Definitions of conditional symmetry for systems of PDEs. 

Here we present new definitions of Q-conditional symmetry which naturally arise for systems 
of PDEs. To avoid possible difficulties that can occur in the case of arbitrary system of PDEs, 
we restrict ourself on the RD systems of the form ([T]). 

It is well-known that to find Lie invariance operators, one needs to consider system ([1]) as 
the manifold A^ = {S'i=0,S'2 = 0} where 

Si = XiUt - u^^ - F{u,v) = 0, 
5*2 = X2Vt - Vxx - G{u, v) = 0, 

in the prolonged space of the variables: t, x, u, v, Ut, Vt,Ux, v^, u^x, Vxx, u^t, v^t, Uu, Vu- According 
to the definition, system ([T]) is invariant under the transformations generated by the infinitesimal 
operator 

Q = i^{t,x,u,v)dt + C.^{t,x,u,v)dx + r]^{t,x,u,v)du + rf{t,x,u,v)dy, (6) 
if the following invariance conditions are satisfied: 

QSi = Q{\iUt - u^x - F{u, v)) 
QS2 = QiX2Vt-v^^-G{u,v)) 

The operator Q is the second prolongation of the operator Q, i.e. 

where the coefficients p and a with relevant subscripts are expressed via the functions r]^ 
and f]^ by well-known formulae (see, e.g., [HI |22l 123]). 



= 0, 

(7) 

= 0. 

M 



9. 



The crucial idea used for introducing the notion of Q-conditional symmetry (non-classical 
symmetry) is to change the manifold Ai, namely: the operator Q is used to reduce Ai. It can 
be noted that there are two essentially different possibilities to realize this idea in the case of 
system ([1]). 

Definition 1. Operator is called the Q-conditional symmetry of the first type for the RD 
system ([T]) if the following invariance conditions are satisfied: 



QSi = Q{XiUt - u^^ - F{u, v)) 
= Q{\2Vt - v^^ - G{u, v)) 



= 0, 
= 0, 

Ml 



where the manifold Aii is either {Si = 0, S2 = 0, Q{u) = 0} or {Si = 0, S2 = 0, Q{v) = 0}. 
Definition 2. Operator ([6]) is called the Q-conditional symmetry of the second type, i.e., 
the standard Q-conditional symmetry (non-classical symmetry) for the RD system ([T]) if the 
following invariance conditions are satisfied: 



^Si = Q{\iUt - u^x - F{u, v)) 
QS2 = Qi\2Vt-v.,.,-G{u,v)) 



= 0, 

M2 ^^Q) 

= 0, 

M2 



where the manifold M2 = {Si = 0,82 = 0, Q{u) = 0, Q{v) = 0}. 

It is easily seen that A^2 C M.i C M., hence, each Lie symmetry is automatically a Q- 
conditional symmetry of the first and second type, while Q-conditional symmetry of the first 
type is one of the second type. From the formal point of view is enough to find all the Q- 
conditional symmetry of the second type. Having the full list of Q-conditional symmetries of 
the second type, one may simply check which of them is Lie symmetry or/and Q- conditional 
symmetry of the first type. 

On the other hand, to construct Q-conditional symmetries of both types for a system of 
PDEs, one needs to solve new nonlinear system, so called system of determining equations, 
which usually is much more complicated than one for searching Lie symmetries. This problem 
arises even in the case of linear single PDE and it was the reason why G.Bluman and J.Cole 
in their pioneering work [21] were unable to describe all Q-conditional symmetries in explicit 
form even for the linear heat equation. Thus, both definition are important from theoretical 
and practical point of view. 

It should be noted that Definition 2 was only used in papers [TOl [251 126] devoted to the 
search Q-conditional symmetries for the systems of PDEs. Moreover, to our best knowledge, 
nobody has noted that a hierarchy of conditional symmetry operators can be defined for systems 
involving, say, m PDEs. In fact, different definitions can be formulated for such systems in 
quite similar way to Definitions 1 and 2 (see [27] for details). 

3. Conditional symmetries of the DLV system 

First of all, we construct the system of determining equations(DEs) to construct Q- 
conditional symmetries of the second type (nonclassical symmetries) of system ([2]). The most 
general form of such operators is the first-order operator 

Q = ^°{t,x,u,v)dt + i^{t,x,u,v)dx + r]^{t,x,u,v)du + r]'^{t, x,u,v)dv, (11) 



where the functions ^'(t, x, u, v) and r]''{t, x, u, v) should be determined from the relevant system 
of DEs . In the case x, u, v) 7^ 0, this system can be reduced to that with ^^(t, x, u,v) = 1 
[27] so that we are looking for the operators 

Q= dt + ^{t,x,u,v)d^ + r]^{t,x,u,v)du + r]'^{t,x,u,v)d^. (12) 

Note we examine system ([2]) only in those case when one is a real system of coupled equa- 
tions, i.e., 62 +^1 7^ 0, and contain non-linear equations (otherwise the system is rather artificial). 

Let us apply Definition 2 to construct the system of DEs for finding operator (fT2!) . According 
to the definition the following invariance conditions must be satisfied: 



Q {XiUt = u^x + u{ai + hiu + Civ)) 

Q {X2Vt = + V{a2 + h2U + C2V)) 



= 0, 

M 

= 0, 

M 



(13) 



where the manifold 

M = {Si= 0, S2 = 0, + ^u, = rj\ vt + Cv. = V^} (14) 
(here the left-hand-sides of ([2]) are denoted as 5*1 and 5*2) while 

^ = Q + Pli; + Pti; + Pl^ + Pli; 

I'^t'ra... I'^t'T'Ki.. i'^ita.... ''^tta... i a., U . 



*^ dutx dvtx ** dutt ** dvtt du^x dvxx 

is the second order prolongation of the operator Q and its coefficients are expressed via the 
functions ^,1]^, and rj"^ by well-known formulae (see, e.g..[TT | [221 123] ). 

Now we apply the rather standard procedure for obtaining system of DEs, using the in- 
variance conditions (fT3|) . From the formal point of view, the procedure is the same as for Lie 
symmetry search, however, four (not two !) different derivatives, say UxxiVxxiUt and Vt-, can be 
excluded using the manifold Ai. After straightforward calculations, one arrives at the nonlinear 
system of DEs 

^uu ^vv ^uv 0, (15) 

Vl = 0, (16) 

Vlu = 0, (17) 

2Aia« + vL - = 0, (18) 

2A2ee. + - 2U = 0, (19) 

(Ai + X2m. + 2r]l - 2^xv = 0, (20) 
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{Xi + X2)CCu + 2vl-2U = 0, (21) 
(Ai - X2)^vl + ^Vlv + 2n(ai + hu + Civ)^, - 2X,C,r]' = 0, (22) 
(As - Xi)^7]l + 2vL + Ma2 + b2U + C2v)^u - 2X2^uV^ = 0, (23) 



Ai(2^„?7i - 6 - ^vV^ - 2^^x) + X2^vV'^ - 3^uu{ai + + av)- 
-ivv{a2 + ^s-u + C2v) - 2r]l^ + S,xx = 0, 

X2{2CvV^ - 6 - CuV^ - 2^Q + Aie^r^i - 3^,via2 + b2U + C2v)- 
-^uu{ai + biu + Civ) - 2r]l^ + ^^x = 0, 

XiiVt + v'^vl + 2^xV^) - X2rf'nl - r?i(ai + 2hiu + civ) - ci'r]2V+ 

+r]lu{ai + biU + Civ) - 2^^-u(ai + biU + Civ) + r]lv{a2 + 62^ + C2v) - r]].^ = 0, 



(24) 
(25) 
(26) 
(27) 



X2{Vt + V'^vl + 2ixrf) - Xir]^r]l - 772(03 + 62^ + 2c2i;) - b2r]iv+ 
+rjlu{ai + biu + civ) - 2^^f (02 + &2'« + C2v) + 'qlv{a2 + b2U + C2v) - rjl^ = 0. 

It turns out, the functions rj^ and 77^ can be at maximum linear functions with respect to 
u and V. In fact, the differential consequences of fl2U]) and fl?I]) with respect to these variables 
lead to the expressions 

(Ai + A2)C' = 0, (Ai + X2)eu = 

so that = = 0. Having ^ = C,(t, x) equations f|T6l) - fl2T|) can be easily solved and one arrives 
at 

= g^(t, x)v + r^(t, x)u + p^{t, x) .^gN 
rf = q^{t,x)u + r'^{t,x)v + p'^{t,x), ^ ' 

hence, the most general form of operator (fT2|) for system ([2]) is as follows 

Q = dt + idx + (q^v + r^u + p^)du + (g^M + r^v + p^)^^, (29) 

where the functions ^ y^^p^ = 2) should be found from other equations. Substituting 



and Cu = = into equations fl221) - fl7n) . those can be splitted with respect to u, v, v?, v"^, uv. 
Finally, one obtains the system of DEs 

(ci - C2)q' = 0, (30) 

(bi - 62)?' = 0, (31) 
Ciq^ + biir' + 2Q=0, (32) 













(33) 


(26i-62)g' + ci(r2 + 2e.) = 0, 










(34) 


(2c2 - c,)q' + b2{r' + 2Q = 0, 










(35) 


(Ai - X2W + 2ql = 0, 










(36) 


(A2 - \iW + 2ql = 0, 










(37) 


Ai(6 + 2eex) + 2r^-e.. = 0, 










(38) 


A2(6 + 2eex) + 2r2-e.. = 0, 










(39) 


Ai(ri + 2rie.) + (Ai-A2)g^g'- 




- 2hp' 




= 0, 


(40) 


A2(r2 + 2r2e,) + (A2-Ai)gig2_ 


h2p'- 


- 2C2P^ 


- 2a2ix - 


= 0, 


(41) 


Ai(g,i + 2gie.) + (Ai-A2)?V- 


(ai - 


' a,2)q^ - 


- ciP^ - qlx = 


0, 


(42) 


A2(gt +2g2ex) + (A2-Ai)gV + 


(ai - 


- a2)q^ - 


■ b2P^ - qlx = 


0, 


(43) 


Ai(p^ + 2piex) + (Ai-A2)(zy- 




Pxx 


0, 




(44) 


A2(p? + 2p2^,) + (A2-Ai)gV- 


a2P^ 


-p' = 

fxx 







(45) 



to find Q-conditional symmetry operator (!29l) of system ([2]). 

Theorem 1 In the case Ai 7^ A2, DLV system ^ is Q- conditionally invariant under operator 
(d^j if and only if 61 = 62 = b, ci = 02 = c. Moreover, z/ foe = 0, 6^ + 7^ then system ^ 
and Q-conditional symmetries of the second type (up to local transformations u — )■ v — >■ 
exp(||t)t>, 6 7^ and u — exp(^t)t', cf — )■ tt, c 7^ ) have the forms 

XiUt = Uxx + u{ai+u), 

X2Vt = fxx + VU, 

and 

Q = dt + T — 9x + (exp(aix + ^t){{a3 + ^4 exp(--^t))it + a-sai) + a2v)dv, (47) 
Ai — A2 V A2 A2 / 



where a^, k = 1,4 are arbitrary constants with the restriction + a1 ^ 0. 

//fee 7^ and the additional restrictions qI = Qx — ^ ^^^^ place then exactly three cases (up 
to local transformations u — )■ hu, v ^ cv and u ^ v, v ^ u ) exist when system ^ admits 
Q-conditional symmetry operators. They are listed as follows: 

/A XiUt = Uxx + u{ai + u + v) , .^g. 
\2Vt = Vxx + v{a2 + u + v), 

Qi = (Ai - X2)dt - {aiv + + 0102) (5„ - d^), aj + aj^ 0, (49) 
Q2 = (Ai - X2)dt + (ai - a2)u{du - d^), ai ^ 02, (50) 

Qs = (Ai - \2)dt - (fli - a2)v{du - d^), ai ^ 03. (51) 

XiUt = Uxx + u{a + u + v), 
X2Vt = Vxx + v{a + u + v), 

Qi = (Ai - X2)dt - a{v + U + a){du -(?„), a 7^ 0, (53) 
Q2 = (Ai - X2)tdt - {Xiv + X2u){du - d^). (54) 

{iii\ XxUt = Uxx^u{aXx^u^v), .^^^ 
X2Vt = Vxx^ 'v[a-X2^ v), a 7^ 0, 

Q\ = (Ai - X2)dt - a{Xiv + X2U + aXiX2){du - dy), (56) 
Q2 = dt + au{du~d^), (57) 
Q3 = dt- av{du - dy), (58) 

Qi = (e""* - a(Ai - X2))dt + aa{XiV + X2U + aAiA2)(9„ - d^), a 7^ 0. (59) 

Proof. Using equations (|5U]) - (PT]) from the system of DEs, one notes that three cases can 
only take place: 1.1 bi 7^ 62 and/or Ci 7^ C2; 1.2 &i = 62 = ^ 7^ 0, Ci = C2 = 0, and 1.3 
61 = 62 = 7^ 0, ci = C2 = c 7^ 0. Note the fourth possible case 61 = 62 = 0, ci = C2 = c 7^ is 
reduced to the second case by the renaming w — )■ f , — )■ u. 

It turns out that case 1.1 produces the restrictions 

g' = 0, q' = 0, (60) 
which lead only to Lie symmetry operators of system ([2]). Let us show this. 
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If bi 7^ &2 and Ci 7^ C2 then restrictions fl60l) immediately follow from fl30l) - fl31l) . If 61 = 62 = ^ 
and Ci 7^ C2 then = follows from (1501) . furthermore, equations and (1531) lead to 
(c2 — ci)q^ = -x^ q'^ = (the subcase 61 7^ 62, ci = C2 = c leads to the same result). 

Having restrictions fl60l) . we immediately obtain cip^ = 0, 62^^ = from fH2|) and fj^3|) . 
Obviously, if C162 7^ 0, then 

pi=j52^0. (61) 

If C162 = 0, say, ci = 0, &2 7^ ( subcase ci 7^ 0, 62 = can be treated in the same way ) 
then = 0. Simultaneously system ([2]) reduced to one with an autonomous equation, which 
is nothing else but the Fisher equation 

XiUt = u^x + u{ai + u), fli 7^ 0. (62) 

Now we assume that such system is Q- conditionally invariant under an operator of the form 
(!29|) . However, the Fisher equation doesn't admit any Q-conditional symmetry [17j but only 
Lie symmetry 

Q = dt + adx, OL = const, (63) 

hence, setting ^ = a into (132|) - (H5|) . we arrive at (I6T1) . The special value ai = in (162|) leads 
to the additional Lie symmetry 2tdt + xdx — 2udu but there aren't Q-conditional symmetry 
operators [T3|. So, the restrictions (161 p are still obtained. 

Restrictions fl6Up and (16 ip simplify essentially the system of DEs (I30p - (1450 . which takes 
the form 

bi{r' + 2Q=0, (64) 

b2{r' + 2Q=0, (65) 

ci(r2 + 2e.) = 0, (66) 

C2(r' + 2ex.) = 0, (67) 

Ai(6 + 2ax) + 2ri.-e.. = 0, (68) 

A2(6 + 2eex) + 2r2-e.. = 0, (69) 

Xi{rl + 2r'^,)-2a,^,-rl = 0, (70) 

A2(r2 + 2r2e.)-2a2e.-rL = 0. (71) 
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Now one may easily check that any solution of system fl64p -( 17T]) leads to a Q-conditional 
symmetry operator which will be equivalent to Lie symmetry operator obtained in [9] . Consider, 
for example, the most general case 62C1 7^ 0. Obviously, fl65|) and fl66|) with 62C1 7^ lead to 

r^ = r^ = -2i^. (72) 

So, having fl72|) and Ai 7^ A2, we obtain the system 

it + 2^^. = 0, 
f =0 

from (l68l) . (|69|) and fl72l) . Substituting the general solution of fl73|) 

e(t,x) = ^ (T4) 

It + a2 



into fl72|) . we can solve equations fl70|) -f l7T|) . Finally, the operator 

Q = dt+ ^'^^ - — {udu + vd^) (75) 
It + 0^2 zt + 0^2 

is obtained if ai = 02 = 0. However, operator (!75l) is nothing else but a linear combination of 
Lie symmetry operators of ([2]) |9] (see table 1, case 1.) multiplied by 2t+a2 ' If ^1 + ^2 7^ then 
operator (163|) occurs, which is Lie symmetry operator. Thus, case 1.1 is completely examined. 

Consider case 1.2. Here system ([2]) can be reduced to one by the substitution u — )■ 
hu, V — )■ exp(||t)t>, because ci = C2 = 0. Since the first equation of ( l46l) is the Fisher equation 
f l62|) the same approach can be used as above. Thus, using operator f l63|) . we again substitute 
i = a into fl32|) - fH5|) . what leads to the restrictions = = = and operator f l29|) takes 
the form 

Q = dt + ad, + {q\ + r\ + p^)d,. (76) 
Simultaneously, the system of DEs fl30|) - fH5l) reduces to 

= = 0, (77) 



(A2 - \iW + = 0, (78) 
X2q^t + - / - = 0, (79) 

A2P?-pL = 0. (80) 
The general solution of this system can be straightforwardly constructed and it reads as follows 

r' = «2, (81) 



P 



g' = c(t)exp(|(Ai-A2)x), (82) 
= exp (|(Ai - A2)a:) (A2c'(t) + aic(t) - (^(Ai - A2))') , (83) 

where 

2 1 ^ 

c(t) = a3exp(^£^(Ai-A2)'t) + a4 exp (^— ( - ai + y (Ai - Aa)')^) , (84) 

and Ofc, A; = 2,3,4 are arbitrary constants. Finally, introducing the notation a = j^zi^'^i 
substituting the function r^jq"^, into fl7B|) . we arrive at the Q-conditional symmetry operator 

dlZD- 

Consider case 1.3. Here the DLV system ([2]) is reduced to system (HHjl . i.e. ([2]) with 
6 = c = 1, by the substitution u — )■ hu, v — )■ cv. Now we take into account the restrictions 
ql = ql = 0, hence, the system of DEs reads as follows 



(Ai - \2)^q' = 0, 










(85) 


(A2 - Ai)^g = 0, 










(86) 


q'+r^ + = 0, 










(87) 


q^ + r' + = 0, 










(88) 


Ai(6 + 2ax) + 2ri- 










(89) 


A2(6 + 2ax) + 2r2- 


^XX 0) 








(90) 




- A2)g^g^ 




2p^ - 2aii^ - 


=0 

' XX 


(91) 




- Xi)q'q' 




2p' - 202^^ - 


^xx 0, 


(92) 


Xi{ql + 2q%) + {X^ 


- \2)q'r' 


- (CLI - 


- a2)q^ -p^ = 


0, 


(93) 


A2(g? + 2g'ex) + (A2 


- Ai)?V 


+ (fll - 


- a2)q^ -p^ = 


0, 


(94) 


X,{pl + 2p'^^) + {X, 


- A2)gy 


- aip^ 


Pxx 0, 




(95) 


\^(p^^+2p^Q + {\2 


- Ai)?V 


- a2P^ 


- pIx = 0. 




(96) 



in 



If = = then we again obtain only Lie's operators (see the case 1.1). So, non-trivial 
results are obtainable only under restriction (q^)"^ + {q^Y ^ 0. Equations under this 

restrictions produce ^ = 0, hence, we obtain 

q^ = = -ip{t), q^ = -r^ = -ip{t) (97) 

from fl87|) - fl90|) (here Lp{t) and '0(t) are arbitrary smooth functions at the moment). Substituting 
(1971) into (I9l-(l9l, we find 



pi = (ai - a2)tP{t) + (A2 - Xi)i^\t) - Xif{t), 
p2 = - aiMt) + (Ai - X2)ipHt) - X^ip'it). 

Having (!97|) and ( |98|) . equations (l9T!) - (p2l) can be rewritten as ODEs for the functions Lp{t) and 

'^'(^) = -T^^T^ («2 - fli + (Ai - \2){v + V')) ((3Ai - Aa)^? + 2A2V') (99) 
1^1 ~ A2J 

^'(^) = 7^^T^(«2 - «i + (Ai - A2)(¥^ + V')) (2Ai(/^ + (3A2 - (100) 

(^Ai — A2J 

Finally, using formulae f pTjl - fllOOl) . the last two equations, fl^ - fl^ . can be rewritten as 
two algebraic equations to find Lp{t) and ip{t). The difference of those leads to the classification 
equation 



(101) 



'fli - 02 - (Ai - A2)(v5 + ^)) + A2V') (ai(4Ai + 5A2) - a2(5Ai + 4A2)- 
-4(Ai - A2)((2Ai + \2)ip + (Ai + 2AA2)^)) = 0. 

Thus, three subcases follow from (IIOII) : 

1.3.1 ^ = -V^ + f3; 

1.3.2 ^ = -Al^; 

1-3-3 V> = 4(Ai-A2K2Ai+A2) ((4^1 - ^"2)Ai + (5ai + 4a2)A2 - 4(Ai - A2)(Ai + 2A2)V 
In subcase 1.3.1, both equations, fl95|) and f l96|) . are equivalent to the equation 

«i - (Ai - A2)^)^(ai - a2 - (Ai - A2)v) = 0. (102) 



If ip = xf^^ then (/? = — xi~Aj ^^"^ using and we arrive at operator 
If ■i/' = 0, then 99 = f^^rfl , so that operator fISUl) is obtained( the restriction ai 7^ 02 guarantees 
that it is no Lie's operator ). 

ffip = f^Efl ; then = 0, what leads to the operator flSTj) and the same restriction oi 7^ 02- 
Thus, the proof of item (i) is completed. 
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In subcase 1.3.2, both equations, (^5^ and OHS]), are equivalent to the equation 



'!/'(t)(ai - a2)(aiA2 - a2Ai) = 0. (103) 

Since ip{t) ^ (otherwise one arrives at the Lie operator Q = dt) two possibihties occur: 
ai = a2 and ai = ^0,2- 

If Oi = 02, then we obtain 

(Ai - A2)t + Aitt 

from f llOOl) . Note we set a = without losing generality. Now the functions p^, p^, can be 
found from (l98i) . hence, we obtain operator (15^ . Operator (1531) follows from (HHl) if one sets 
ai = 02 = a. Thus, the proof of item (ii) is completed. 

If ai = j^a2 (here 02 7^ otherwise item (ii) is obtained) then equation (llOOp produces 

q;o02Ai 



exp(-f|t) - ao(Ai - A2)A2 



where is a non-vanish constant. Thus, using equations (p8|) and notations 02 = C1A2, A2ao = 
a, we obtain the most complicated operator (15^ . Finally, operators (15^ . and are 
nothing else but those fH^ . (15(1 and (IFI]) with Oi = aAi,a2 = aA2, respectively. Thus, all 
operators arising in item (Hi) are constructed. 

It turns out that the detailed analysis of subcase 1.3.3 doesn't lead to any new operators. 
The proof is now completed. ■ 
Remark 1. If the restrictions ql = = don't take place we were not able to solve 
the corresponding nonlinear system of DEs, hence, DLV system ([2]) may admit Q-conditional 
symmetries of the form fl2I?l) with ql and/or 7^ 0. 

Theorem 2 In the case Ai = A2, DLV system ^ admits only such operators of the form ^2^, 
which are equivalent to the Lie symmetry operators. 

Theorem 3 In the case Ai 7^ A2, DLV system ^ is invariant under Q-conditional operators 
of the first type only in two cases. The corresponding systems and Q-conditional symmetries 
(up to local transformations u bu, f — )■ exp{^t)v, b ^ and u — )■ exp(|^t)t;, cv u, c ^ 0) 
have the forms 

/A XiUt = u^:r + u{ai + u + v), 

Qi = (Ai - X2)dt + (ai - a2)u{du - dy), (105) 
Q2 = (Ai - X2)dt - (tti - a2)v{du - dv)- (106) 

/-A XxUt = Uxx^u{ax^u), 
X^vt = + VU, 
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Q = dt + - — ^^^dx+ ( exp{aix + ^t) [{as + a4exp{-^t))u + a3ai) +a2v)dy, (108) 
Ai — A2 V A2 A2 ^ 

where ak, A; = 1, ... ,4 are arbitrary constants with the restriction al + al ^ 0. There are no 
any other Q-conditional operators of the first type. 

In the case Ai = A2, DLV system ^ is invariant only under such Q-conditional operators 
of the first type, which coincide with the Lie symmetry operators. 

Proofs of Theorems 2 and 3 are similar to one presented above for Theorem 1 and omitted 
here because their bulk. It should be noted that both manifolds arising in Definition 1 and the 
most general form ([TT]) of the operator in question were used. 

Remark 2. Theorems 2 and 3 give a complete description of Q-conditional symmetries of 
the first type in explicit form because there aren't any additional restrictions on the form of 
those operators (in contrary to the Q-conditional symmetries of the second type). 



4. Reductions to ODEs' systems, exact solutions and their application 



First of all, we note that DLV system ([2]) is invariant under time and space translations, 
hence, its arbitrary solution uo(t,x), Vo{t,x) generates a two-parameter family of solutions of 
the form uo(t — tQ,x — Xq), Vo{t — to,x — Xq). Having this in mind, we set to = = in the 
solutions obtained below. 

It is well-known that using Q-conditional symmetries one can reduce the given two- 
dimensional PDE (system of PDEs) to an ODE (system of ODEs) via the same procedure 
as for classical Lie symmetries. Thus, to construct an ansatz corresponding to the operator Q, 
the system of the linear first-order PDEs 

should be solved. Substituting the ansatz obtained into DLV system with correctly-specified 
coefficients, one obtains an system of ODEs, i.e., the reduced system of equations. Since this 
procedure is the same for all operators, we consider only operator (j^9l) in details. In this case 
system (11091) takes the form 

{Xi - X2)ut = -{aiv + a2U + aia2), 
{Xi - X2)vt = aiv + + aia2. 

To solve flllOp we immediately note that Ut = —Vt, hence, 

u{t,x) = —v{t,x) + ipi{x). (Ill) 
Substituting (lllip into the second equation of (1 11 01) . we arrive at the equation 

(Ai - X2)vt = (ai - a2)v + a2V?i(x) + aia2. 
If ai 7^ 02, then this equation has the general solution 

v{t, x) = — ^ ( exp( -^-i^ — '^t)ip2{x) - a2(pi{x) - 0102) , 

(3,1 — a2 Ai — A2 

therefore the ansatz 

^(^'^) = ^73^( -exp(f^t)(^2(a^) + ai¥'i(a;) + 0102), 
^(^'^) = ^7^(exp(f^t)¥?2(2;) -a2V5i(a;) -aitta) 

is obtained. Here (fi and Lp2 are to be found functions. 
If Oi = 02 = a, then the ansatz 

u{t, x) = ipi{x) - if2{x) - j^ifiix) + a)t, 
v{t, x) = (p2{x) + j^iMx) + a)t 

is obtained. 

U 



(112) 



(113) 



To construct the reduced system, we substitute ansatz flll2p into 0481) . It means that we 
simply calculate the derivatives Ut, Vt, u^x, v^x, and insert them into (HSj) . After the relevant 
simplifications one arrives at the ODEs system 

ip'l + (pl + (ai + a2)(pi + aia2 = 0, 

(114) 

,J' _|_ a2Ai-aiA2 ,„ _i_ ,n ,n D 

^2 A1-A2 + ^1^2 — U 

to find the functions (fi and (/?2- Similarly, ansatz flllSp leads to the reduced system of equations 

ifi'l + ia + ipif = 0, 

v'2 + i^2-^^){a + ^,) = 0. ^^^^^ 

In a quite similar way other operators listed in Theorem 1 were used to find ansatze and 
reduced systems of ODEs. They are presented in Table 1. 

Now we construct exact solutions of DLV system using the ansatze and the reduced systems 
obtained above. It should be stressed that all the ODE systems listed in Table 1 are nonlinear 
and none of them can be easily integrated. 

Let us consider system fllMp obtained by application of ansatz flll2p . Since the general so- 
lution of this nonlinear ODE systems cannot be found in an explicit form, we look for particular 
solutions. Setting Lpi = a = const, we find 

+ (oi + a2)a + 0102 = ^ ai = — ai, a2 = —02 

from the first equation of system (11140 . Now we take Lpi = — Oi (the case y?i = —02 leads to 
the solution with the same structure) and substitute into the second equation of system (I114p : 

- (3X,ip2 = 0, (116) 

where = f^^rf| 7^ 0. Depending on sign of the parameter /3 the linear ODE (I116p generates 
two families of the general solutions. Using those solutions and ansatz (11121) . we obtain the 
following two families of exact solutions of the DLV system 



u{t, x) = -ai + (Ci exp{y/^x) + C2 exp(-y^a;)) 
vit, x) = ^ (Ci eM^fWix) + C2 eM-^fWix))e^\ ^^^^^ 

if /3 > 0, and 

u{t,x) = — fli H ^(ClCos{^/^'^X^x) + C2 sin(v^— /3Aiz))e''*, 

a2 ai \ J (1181 

(t,x) = ^(Cicos(v/q5A7x) +C2sin(v/=;5A7x))e^*, 



V 



if /? < (hereafter Ci, C2 are arbitrary constants). 



Table 1. Ansatze and reduced systems of ODEs for DLV system ([2]) 







Ansatze 


Systems of ODEs 


1 


'1 !' 


\ ' J rJ-V/" J- 

x) = 992(w)e^2t + (^^^^Cio; + At) x 
((Cs + C4exp(fit))v.iH + aiC4exp(fii)) 


ip'^ + Ci Ai ip\ + (Oi + ip^ )ip-i = 
(/72 + CiX2<f'2 + ^2i<fl - C2X2) = 


2 


(m 


^(i,x) = ^(-exp(«tV2(x)+ 
a\^\{x) + 0102) 

a2(^i(x) - 0102) 


(/?'/ + V'f + (ai + 02)991 + 0102 = 

,„// 1 12^1 — aiA2,„„ 1 n 

^^2 + A1-A2 ^2 + ¥'l'/'2 — U 


o 




u{z,x) — (p2{x)eW\\-^-\^i^) 
v{t,x) = (pi{x) (p2{x)exp{^]_^lt) 


/o" — U ^/n^ — L rtr^in-t — 

^2 + "1::r/V2 + v^iy^2 = o 


4 


m 


u{t,x) = ipi{x) - ^2ix) expi^^^t) 
v(t,x) = ip2ix)exp{^t) 


(/i?l +1/21+ Oi<y9i = 
^2 + "1::r/V2+ 9^1^2 = 


5 




u{t,x) = ipi{x) ip2ix) ^^^^^(ipiix) + a)t 
v{t,x) = (p2{x) + ^^'^^^(ipiix) + a)t 


v?;' + (o + (^i)2 = o 

'^2 + (<^2- Af^)(o + ^l) = 


6 


(IMl) 


u{t,x) = ^ i;^^ {\wi{x) tLp2{x)) 

■^(*'^) = A,-A, ( >y2Vl{x) + tip2{x)) 


(/^'l' + (/52 + (/7i (0 + y?i) = 

(/?2 + ip2{a + (^1) = 


7 


([56]) 


^(*'^) = a(Ai-A2)( e"V2(x) + 

aXufiix) — a?\i\2) 

^(t,x) = aix'-X,)i^"'Mx) 

aX2(fi{x) - 0^X1X2) 


+ + o(Ai + X2)ip\ + o2AiA2 = 

(/?2 + </?lV?2 = 


8 


dSZl) 


u{t, x) = e"-^(p2{x) 
v{t,x) = ipi{x) — e°'*(p2{x) 


¥'1 +</'!+ «A2(/?1 = 
(/7'2' + V'lV'2 = 


9 




u{t,x) = (pi{x) - e°'^(f2{x) 
v{t,x) = e"-^ip2{x) 


y?" + + oAiv?i = 
993 + (^i9J2 = 


10 




u(t,x) = j^^{Xiipi{x) + aAiA2- 
V?2(x)(l-a(Ai-A2)e'^*)) 
v{t, x) = ^(v,2(x)(l - a(Ai - A2)e»*)- 
X2(pi{x) - aXiX2) 


'/'l + V^l - 0^52 + 

o(Ai + A2)<^i + o2AiA2 = 

V92 + V3i(/?2 = 



Remark 3. In Table 1, the parameter A = -^^Cf — while C^, k = 1, ... ,4 are 
arbitrary constants. 



Let construct solutions of f lll4p with some restrictions on Ai and Ai. Firstly, we note that 
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the substitution 

V?i = v5-ai (119) 
simplifies the first equation of flll4p to the form 

^" + + (as - ai)(^ = 0. (120) 
Of course, (11201) can be reduced to the first-order ODE 

= -l^' + {a,-a,)^' + C (121) 

with the general solution containing special functions, Weierstrass functions [28]. To avoid 
cumbersome formulae, we set C = 0, hence, the general solution is 

^ = 2^0'! - '^2)(l - tanh^(-Vai - 0.2 x)), (122) 

3 I 
if = -(ai - 02) (1 - coth^(-A/ai - a2x)), (123) 

if tti > a2, and 

3 2 

ip= -{ai- a2)(l + tan^(-A/a2 - ai x)), (124) 
if tti < a2. 

Thus, we can apply each of formulae fll22p - fll24l) to solve the second equation of f lll4p . In 
the case of solution f ll22p . this ODE takes the form 

cp2 + (/?2(ai - aa) {^jj^ZJ^ ~ ^ tanh^ {^^/a^^^ x)) = 0. (125) 

Nevertheless, the general solution of fll25p is still unknown, its particular solutions can be found 
[29]: 

„ 1 

ip2 = cosh (-V^i - a2x), (126) 



if Ai = |A2, and 



ip2 = sinh(--\/ai — 02 x) cosh'^ ( - V^i — a2x), (127) 
if Ai = IA2. 

Thus, substituting the functions fi{x) and '^2{x) given by formulae (I119p . (I122p and (11260 
into ansatz (I112p . we find the exact solution of DLV system (HHjl 

u{t,x) = ^ai — |ai tanh^dv^ai — 02 x) — 

cosh^d Vai - a2 x) exp(^^^^t), ^^^^^ 

v{t, x) = -§02(1 - ta■D.h^{^^/al - a2x)) + 

+ ^ cosh=^(|v/^r^x) eM'-^t). 
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if Ai = |A2, ai > a2. 

Similarly, solution fll27p leads to the exact solution 



u{t,x) = ^fli — |ai tanh^(|^ai — 02 x) — 

"^73^ sinh(iA/ai — 022;) cosh'^ {^^/al^^a^ x) exp( ^*-"^^~"^'' t), 



x) = -|a2(l - tanh (|a/'^i ~ o-2x)) + 

+ smh{^y/ai — 02 x) cosh^ ( | -\/ai — 02 x) exp( '^^"^^~"^-' t) 



(129) 



of DLV system (148|) with Ai = |A2, ai > 02. 

In a quite similar way solutions f ll23p and f ll24p were also used to construct three new 
solutions of DLV system fH5|) . Omitting straightforward calculations we present only the result: 



u{t,x) = !«! — |ai coth^(|^ai — 02 x) — 

-^73^ sinh^dVai -022:) exp(5^^^t), 
x) = -§02(1 - coth^(|A/ai - 02 

/t; 73 ovT^i'M^'i3^*2) 1 



(130) 



^^1^ sinh^(ly^r^x) exp(^i^t) 



if Ai = IA2, ai > 02; 



5 

u{t,x) = |ai — |ai coth^(|Y^ai — 02 x) — 

-^73^ cosh(iV«i - ^22;) sinh^dV^i - ^2 a;) exp( ^^"\~°'^ t), 
z;(t, x) = -§02(1 - coth^(|v^ai - 02 x))+ 

+ ^^3^ cosh(iY^ai — a2 x) sinh^(|Y^ai — 02 x) exp( '^*-"^^~"^-' t), 



:i3i^ 



if Ai — |A2, cii > 0,2] 



u{t, x) = |ai + |ai tan^(|A/a2 — ai x) — 

cos^dV^^x) exp(^to^t), 
v{t, x) = -§02(1 + tan^(|A/a2 - Oi x)) + 

ai_a2 cos3(iVa2 - ai x) exp(^^^^t), 



(132) 



if Ai = |A2, ai < 02. 

In a similar way one may use other ansatze and reduced systems of ODEs for constructing 
exact solutions of DLV system (|2]) with the correctly-specified coefficients. Let us consider the 
most cumbersome ansatz and reduced system, which correspond to the Q-conditional operator 
( 159|) . Nevertheless the reduced system of ODEs (see case 10 of Table 1) is again non-integrable, 
its particular solutions can be derived by setting 

ip2 = Aiv?i + aAiA2. (133) 
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The reduced system under condition (11331) takes the form 

+ + aX2ipi = 0. (134) 

Since ODE fll34p has the same structure as (11201) . we immediately obtain its solutions (I122p - 
(11241) with a2 — ai = a\2. Thus, substituting (I133P and (I122p - (ll24p with 02 — fli = (1X2 into 
the ansatz listed in the last case of Table 1, we find the exact solutions of DLV system (1551) 

u{t^ x) = aaXiX2{ — | + | tanh^(i^/— aA2 a:))e'^*, 

v{t, x) = -faAa (l - tanh^(| V-aAa x)) - (135) 

-Q;aAiA2(-| + ftanh^ (|V-aA2 x))e''\ 
u{t, x) = aa\i\2{ — ^ + I cotl?{\\/~a\2 x))e"*, 

v{t,x) = -|aA2(l - coth^ {^y/^aX^x)) - ^^35) 
-aaAiA2( - | + | coth^(i V-aA2 a;))e''*, 

if a < and 

u(t,x) = — aaAiA2(^ + ^t&n'^(^^/aX^x))e"■\ 

(137) 

v{t, x) = — |aA2(l + tan^(|A/aA2 x)) + aaAiA2(| + | tan^(iv^aA^x))e°*, 
if a > 0. 

It should be noted that all the solutions obtained above cannot be constructed using Lie 
symmetries. In fact DLV systems ( HHl) and ( |55l) don't admit any non-trivial Lie symmetry, 
hence, the standard plane wave solutions only can be found by Lie symmetry reductions. 

Finally, we present an example that demonstrates remarkable properties of some solutions 
presented above. 

Example. Consider solution (IllSp with Ci = 0. Using substitution u = —bU, v = 
—cV {b > 0, c > 0), one transforms DLV system (HHl) to the system describing the competition 
of two species 

X-^Ut = U^^ + U{ai -bU - cV), 
X2Vt = + V{a2 -bU- cV) 

and solution (IllSp to the form 



(138) 



V{t,x) = C2^m{^^,x)eP\ 



(139) 



where the coefficient restrictions [3 = j^z^ < 0, ai > 0, 02 > are assumed. Using this 
solution one may formulate the following theorem giving the classical solution for the BVP 
with the constant Dirichlet conditions on the boundaries. 

Theorem 4 The classical solution of boundary-value problem for the competition system I^138i) 
and the initial profile 

= f+ C2sin(y^x), 
^(0'^) = C2sin(v/^x), 

IP 



X 



Figure 1: Solution f lT39|) of system ffT38|) with = 1, 02 = 2, Ai = 11, A2 = 1, b = 0.1, c = 
0.1, C2 = 0.2, (3 = -0.1. 

and boundary conditions 

X = (} ■ U = ^ V = 

in domain Q = {{t,x) E (0, +00) x ^0, --^^^=J } is given by formulae ^3E)■ 
The solution f ll39p with /3 < has the time asymptotic 

{U,V)^{^,0), (142) 

Thus, this solution describes the competition between the two species when the species U 
eventually dominate while the species V die. An example of this competition with the correctly- 
specified coefficients is presented on Fig.l. 

5. Conclusions 

In this paper Q-conditional symmetries of the diffusive Lotka-Volterra system ([2]) and their 
application for finding exact solutions are studied. Following the recent paper [27] , two different 
definitions of such symmetries are used to derive systems of DEs and to construct them in the 
explicit form. It turns out that Q-conditional operators of the first type can be derived much 
easier than those of the second type (nonclassical symmetries), hence. Theorem 3 was proved, 
which completes description of Q-conditional operators of the first type. Nevertheless Theorem 
1 presents a wider list of Q-conditional operators (because all the operators of the first type 
are automatically among those of the second type), the result is still incomplete. In fact, the 
additional restriction on the form of operators in question was used. 

All the Q-conditional operators obtained were used to construct non-Lie ansatze and to 
reduce DLV system ([2]) to the corresponding systems of ODEs. As result, a wide range of new 
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exact solutions in explicit form was found. These solutions possess a complicated structure and 
cannot be found by the classical Lie algorithm. In the particular case, they differs from the 
standard plane wave solutions obtained in j9l[30]. Note plane wave solutions can be derived from 
some ansatze listed in Table 1 under additional restrictions (see, e.g., ansatze with (fi = const 
in 3-rd and 4-th cases of the table). 

Finally, a realistic interpretation for competing species has been provided for exact solution 
f ll39p . In fact, this solution describes the competition between two populations of species when 
one of them eventually leads to the extinction of other. 

The work is in progress to construct Q- conditional symmetries and exact solutions of a 
multicomponent DLV system. 
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